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Abstract 

■ 

In this brief note some comments about the results presented in a 
recently published paper on the measurement of the general relativistic 
Lense-Thirring in the gravitational field of the Earth are presented. It 
turns out that, among other things, the authors might have yielded an 
optimistic evaluation of the error budget because of an underestima- 
tion of the impact of the secular variations of the even zonal harmonics 
of the geopotential. More tests with real data by varying the obser- 
O vational time spans and the magnitudes of J4 and Jq in the orbital 

cr 

processor's force models should have been performed in order to cor- 
Vh^, rectly address this important issue. Preliminary analytical evaluations 

point towards a 15-45% error at 1-3(7 level, respectively. 



1 On the adopted observable 



In a recently published paper 1, submitted 2 June 2004 and accepted 10 
September 2004, a measurement of the general relativistic Lense-Thirring 
effect in the gravitational field of the Earth with the laser-ranged satellites 
LAGEOS and LAGEOS II is reported. The claimed total accuracy is 5 — 
10%. The observable used in the analysis is the following linear combination 
of the residuals of the longitudes of the ascending nodes O of LAGEOS and 
LAGEOS II 

^LAGEOS + Cl ^LAGEOS II „ (1) 

in which c\ = 0.546 and //lt is a scaling parameter which is 1 in the Ein- 
stein's General Theory of Relativity and in the Newtonian mechanics. Eq. 
(0 allows to cancel out all the static and time-dependent contributions of 
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the first even zonal harmonic coefficient J2 of the multipolar expansion of the 
Earth gravitational potential which represents one of the major sources of 
systematic error. The terrestrial gravity model adopted in pQ| is the recently 
released GRACE-only EIGEN-GRACE02S model 0. 

The possibility of using only the nodes of the LAGEOS satellites in view 
of the improvements in the Earth gravity field solutions from the CHAMP 
and GRACE missions was first presented in although not in an analytic 
and explicit form. Eq. (^Q) was explicitly published for the first time in [I]. 
Then, it has been discussed in a number of other papers (SHOIIZI- 

Instead, eq. (0) is presented in as an own result of the authors who, 
not only miss out to correctly cite the appropriate works IU El U\ , but 
improperly refer to [SJ (reference 19, pag. 959 of ). In that paper, which 
is almost twenty years old, there is no mention of eq. (0). It is devoted to 
the well known LAGEOS III mission in which the launch of a LAGEOS-type 
satellite in an orbit with supplementary inclination with respect to LAGEOS 
is presented. The goal of that configuration was to cancel out exactly the 
contributions of all the even zonal harmonics of the geopotential by using 
the simple sum of the nodes as observable. 

2 Some possible criticisms on the error budget 

In this Section we will show that the error analysis presented in pQ could 
be considered, perhaps, too optimistic, mainly with respect to the impact of 
the gravitational errors which represent the major source of systematic bias 
when eq. (^Q) is adopted. 

In Section Error assesment, pag. 960 of the authors correctly as- 
sess the systematic error due to the static part of uncancelled even zonal 
harmonics jf^ of geopotential: indeed, it is 3—4%, according to the EIGEN- 
GRACE02S model. This result agrees with that obtained in It should be 
pointed out that these evaluations are at 1 — a level; at, say, 3 — a we would 
get 9 — 12%. The first number (3%) comes from a root-sum-square calcula- 
tion while the second number (4%) is the upper bound obtained by simply 
summing up the individual error terms. The calibrated standard deviations 
of the even zonal coefficients of EIGEN-GRACE02S have been used. A pos- 
sible criticism is that only this Earth gravity model has been used in the pre- 
sented analysis. Using different GRACE-only Earth gravity models, like, e.g. 
EIGEN-GRACE01S (available at http://op.gfz-potsdam.de/grace/results) 
and GGM01S jlUj . and analyzing the scatter of the so-obtained results would 
have yielded better results in term of confidence and reliability. Moreover, 
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also the risk of using that particular model which gives just the expected 
result would have been avoided. 

Other problems may arise when the authors show their a priori error 
analysis for the time-dependent gravitational perturbations (solar and lunar 
Earth tides, secular trends in the even zonal harmonics of the Earth's field 
and other periodic variations in the Earth's harmonics). Indeed, they claim 
that, over an observational time span of 11 years, their impact would be 2%. 
This evaluation is based on reference 30, pag. 960 of which refers to the 
WEBER-SAT/LARES INFN study; it has nothing to do with the present 
node-only LAGEOS-LAGEOS II combination. On the contrary, many re- 
cent studies ^3 El El IH1 El Ej mainly focussed on the gravitational 
part of the error budget in the performed or proposed Lense-Thirring tests 
with LAGEOS-like satellites are not even included in the attached .doc file 
which should overcome the unavoidable space limitations posed by the Let- 
ter format 1 . Moreover, this estimate may turn out to be optimistic because 
of the secular variations of the even zonal harmonics 2 Jg. Indeed, eq. (J2) 
allows to cancel out J2, but is sensitive to J4, Jq,..., as pointed out in [S]. 
The uncertainties in the Jg are still quite large. On the other hand, their 
impact on the Lense-Thirring measurement grows linearly in time. Indeed, 
the mismodelled shift, in mas, of eq. Q due to the secular variations of the 
uncancelled even zonal harmonics can be written as 

£ (r> L / GEOS + Cl r> L / GEOS ") ^ri s , (2) 

where the coefficients Clj are <90 c i ass /9J£ and have been explicitly calculated 
up to degree £ = 20 in JH]- It must be divided by the gravitomagnetic shift, 
in mas, of eq. over the same observational time span 

^LAGEOS + M6 ^LAGEOS ifj ^ = ^ mas ^ ^ (3) 

By assuming c) J4 = 0.6 x 10 -11 yr -1 and 5Jq = 0.5 x 10~ n yr _1 [IS], it 
turns out that the percent error on the combination eq. Q grows linearly 
with T Q bs and would amount to 1% over one year at 1 — a level. This means 

On the contrary, a large number of references are devoted to the non-gravitational 
perturbations which, instead, play a minor role in this case due to the small senistivity of 
the LAGEOS nodes to them. 

2 The problem of the secular variations of the even zonal harmonics in post-Newtonian 
tests of gravity with LAGEOS satellites has been quantitatively addressed for the first 
time in |17|. In regard to the Lense-Thirring measurement with eq. (Q, it has been, 
perhaps, misunderstood in 
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that, over 11 years, their impact might range from 11% (1-cr) to 33% (3-<r). 
Alternatively, if we look at the rate 3 , in mas yr _1 , these figures must be 
doubled. Indeed, the mismodelled secular rate due to the Ji is 

£ (f2 L / GEOS + Cl r2 L / GEOS n ) 5J e T ohs , (4) 

1=2 

which must be divided by the Lense-Thirring secular trend 

rjLAGEOS + M6 ^LAGEOS II = ^ mag yr ~l (5) 

This subtle and important point should have been addressed with tests with 
real data by varying the magnitudes of J% and J4 in the force models of the 
orbital processor over different observational time spans. 

Another controversial point is that it is unlikely that the various errors 
of gravitational origin can be summed in a root-sum-square way because 
of the unavoidable correlations between the various phenomena of gravita- 
tional origin. It would be more conservative to add them. In this case, the 
{jf ] - j e ) error would range from 15% (4%+ll%) at 1 - a level to 45% 
(12%+33%) at 3 — a level over a 11-years long observational time span 4 . 
The so obtained global gravitational error can be added in quadrature to 
the non- gravitational error. Even by assuming the 2% authors' estimate of 
the time-dependent part of the gravitational error, the upper bound errors 
would be ■y/(4 + 2) 2 + 2 2 % = 6% at ^(8 + 4)2 + 42% = 13% at 2-a 
and y 7 (12 + 6) 2 + 6 2 % = 19% at 3-cr. Instead, at the end of the Section 
Total uncertainty, pag. 960 of and in their Supplementary Information 
.doc file the authors add in quadrature the doubled error due to the static 
part of the geopotential (the 2 x 4% value obtained from the sum of the in- 
dividual error terms), their perhaps optimistic evaluation of the error due 
to the time dependent part of the geopotential and the non-gravitational 
error getting \/8 2 + 4 2 + 4 2 % = 10% at 2-a. On the other hand, in the 
Supplementary Information .doc file it seems that they triple the 3% error 
due to the static part of the geopotential obtained with a root-sum-square 
calculation and add it in quadrature to the other (not tripled) errors getting 
V9 2 + 2 2 -|- 2 2 % < 10% at 3-cr. These calculations look like tricks to get just 
a desired value, i.e. 10%. 

Finally, it is hard to understand why the authors very often refer to 
the LAGEOS-LARES proposed experiment and to the related simulations 

3 Indeed, the normalized slope of the time series is measured. 

4 The evaluations of eq. © have been used; if eq. (|1J is used the upper bounds become 
4% + 22% = 26% (1-a) and 12% + 66% = 78% (3-cr). 



4 



Table 1: Orbital parameters of the existing LAGEOS and LAGEOS II and 
of the proposed LARES. 



Satellite 


a semimajor axis (km) 


eccentricity e 


inclination i (deg) 


LAGEOS 


12270 


0.0045 


110 


LAGEOS II 


12163 


0.0135 


52.64 


LARES 


12270 


0.04 


70.0 



and error budgets. It is rather confusing and misleading. The LAGEOS- 
LAGEOS II combination of eq. (JTJ) is, by construction, designed in order 
to exactly cancel out the J% term with an approach which can be applied 
to any orbital configuration given a pair of satellites in different orbits or 
a pair of different Keplerian orbital elements of the same satellite. On the 
contrary, the observable originally proposed for the LAGEOS-LARES mis- 
sion is the simple sum of their nodes. If the orbital parameters of LARES, 
quoted in Tabled were exactly equal to their nominal values, all the even 
zonal harmonics would be exactly cancelled out. Instead, the sum of the 
nodes would be affected, to a certain extent, by the whole range of the even 
zonal harmonics of the geopotential due to unavoidable departures from the 
LARES nominal configuration because of the orbital injection errors and 
mission design (the eccentricity of LARES would be one order of magnitude 
larger than that of LAGEOS), i.e. the coefficients of the classical nodal pre- 
cessions would not be exactly equal and opposite f^ AGEOS ^ — f^ ARES for 
I = 2, 4, 6, 8, .... In the Supplementary Information .doc file the combination 
of the nodes of LAGEOS and LAGEOS II of eq. (0) is presented as if it is 
only slightly different with respect to the sum of the nodes of the originally 
proposed LAGEOS-LARES configuration, apart from a 18 deg offset in the 
inclination of LAGEOS II with respect to LARES. The differences in the 
eccentricities and the semimajor axes, which do play a role [T5|, have been 
neglected. 

3 Conclusions 

The main objections to the results presented in [J can be summarized as 
follows 

• The authors attribute to themselves the combination of eq. Q used 
in their analysis and consistently ignore almost all the works of other 
scientists on the gravitational part of the error budget which, in this 
case, represents the major source of systematic bias because of the 
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low sensitivity of the nodes of the LAGEOS satellites to the non- 
gravitational perturbations. 

• No different GRACE-only Earth gravity models have been used to 
support the error assessment. 

• A major point is represented by the impact of the secular variations 
of the uncancelled even zonal harmonics J4, Jq. Their mismodelling 
might induce errors as large as tens percent over an observational time 
span of 11 years, according to preliminary analytic analyses. Tests 
with real data by varying the observational time spans and the mag- 
nitudes of J4, Jq in the orbital processor's force models should have 
been extensively conducted. The authors present a 2% estimate of the 
global time-dependent part of the gravitational error and support it 
with an improper reference to the preparatory study for the WEBER- 
SAT/LARES mission submitted to the Italian National Institute of 
Nuclear Physics (INFN). 

• The static and time-dependent gravitational errors should be simply 
added together because of the correlations among the various gravi- 
tational phenomena. Instead, in [Q they are summed in quadrature 
and presented in a way which seems a ad hoc trick just to get a 10%. 
Even by assuming the authors' 2% estimate of the time-dependent part 
of the gravitational error, the 3-<r upper bound error would amount 
to 19%, contrary to the 10% value presented in the Supplementary 
Information .doc file. 
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